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A procedure  is  outlined  for  analyzing  shells  of  revolution  with  nonaxisym- 
metric  geometrical  or  physical  occurrences,  e.g.  cutouts.  Such  a cutout  means 
that  the  shell  is  nc  longer  a shell  of  revolution.  Advantage  is  taken  of  the 
efficient  tools  available  for  analyzing  shells  of  revolution  to  find  the  res- 
ponse of  a shell  which  can  be  modeled  as  a shell  of  revolution  with  such  non- 
axlsymmetric  occurrences  as  cutouts  or  attachments.  In  particular,  the  proposec 
scheme  is  a postprocessor  for  the  numerical  integration  method  of  shell  of 
revolution  analysis.*:  A L4 


— 


INTRODUCTION 


A procedure  for  analyzing  shells  of  revolution  with  nonaxisymmetric 


geometrical  or  physical  occurrences,  e.g.  cutouts,  is  outlined.  This 


means,  of  course,  that  the  shell  is  not  actually  a shell  of  revolution 


The  goal  is  to  take  advantage  of  the  efficient  tools  available  for 


analyzing  sheila  of  revolution  to  find  the  response  of  a shell  which  can 


be  modeled  as  a shell  of  revolution  with  such  nonaxisymmetric  occurrences 


as  cutouts  or  attachments 


In  some  stress  analysis  problems,  certain  variables  can  be  expanded 


as  Fourier  series  permitting  the  unknown  Fourier  coefficients  to  be  determined 


such  that  the  governing  differential  equations  and  boundary  conditions 


are  satisfied.  Thus,  in  a three-dimensional  problem  the  introduction  of 


a Fourier  series  into  the  governing  differential  equations  can  eliminate 


Thereby,  a three-dimensional  problem 


one  coordinate  from  the  probl 


can  be  reduced  to  a two-dimensional  problem.  It  is  in  this  manner  that 


the  circumferential  coordinate  is  frequently  eliminated  from  the  shell  of 


revolution  governing  differential  equations.  The  resulting  set  of  simultaneous 


ordinary  differential  equations,  rather  than  the  original  partial  differential 


equations,  are  often  solved  using  finite  differences  (e.g.  (1))  or  the 
nine r leal  integration  (transfer  matrix)  method  (e.g.  (2,3,4)  and  numerous 
subsequent  papers  by  the  same  authors).  However,  if  there  is  a nonsymmetric 


cutout  or  attachment  this  method  can  not  be  applied  directly.  In  this 


technical  note,  a combination  of  the  mnerlcal  Integration  method  and  the 


finite  difference  method  could  be  used  rather  than  the  finite  element 


method  to  couple  to  the  numerical  Integration  technique. 

For  the  scheme  developed  here,  the  numerical  Integration  method  Is 
proposed  for  region  I and  III  of  Fig.  1 and  the  finite  element  method  can 
be  used  for  region  II.  To  employ  the  combination  of  the  two  methods,  we 
need  to  transform  the  field  (transfer)  matrix  developed  by  numerical 
Integration  into  a stiffness  matrix  or  the  stiffness  matrix  of  the  finite 
element  method  into  a transfer  matrix.  If  the  nodal  variables  in  the 
transfer  matrix  and  the  stiffness  matrix  are  of  the  same  type  as  each 
other,  this  can  be  accomplished  using  a simple  rearrangement  of  variables 
(5).  However,  since  for  the  numerical  integration  the  variables  are 
expanded  as  a Fourier  series,  this  method  produces  a state  vector 
containing  Fourier  coefficients  rather  than  the  physical  variables  (e.g., 
displacements  themselves).  Therefore  it  is  necessary  to  make  the  variables 
compatible  with  the  nodal  physical  variables  of  the  finite  element 
method.  This  can  be  accomplished  with  the  following  formulations. 

FORMULATION 

Let  _z  be  a state  vector  which  includes  generalized  displacements  and 
generalized  forces  at  a station.  Also,  let  x be  the  longitudinal 
coordinate  and  9 be  the  circumferential  coordinate  of  a shell  of 
revolution.  Then 

z - z(x,6)  (1) 

If  state  vector  £(x,0)  is  expanded  as  a Fourier  series  (6) 

00 

z(x,9)  ■ E [zc(x)cosmO  + zs(x)slnm0]  (2) 

m-0  " ■ 

£ 

in  which  z (X)  is  a vector  of  coeiiicients  of  a Fourier  cosine  series  and 
z*(x)is  a vector  of  coefficients  of  a Fourier  sine  series.  For  the  purpose 

of  demonstrating  the  technique,  consider  only  m-0,  1 and  0-0,  0^,  The 
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formulation  is  readily  extended  to  more  general  cases,  i.e.  more  terms 
in  the  series  and  more  circumferential  legations. 


Derivation  of  the  stiffness  matrix  from  the  transfer  matrix.  - Let  x^  » 0 
and  x^  correspond  to  the  longitudinal  coordinates  of  the  edge  of  region  I as 
shown  in  Fig.  1.  If  x - Xq  and  x ■ x^  are  introduced  into  Eq.  2, 

i:(x0,e)  - £ [z^(xQ)  cosm0  + z^(x0)  sinm0]  (3a) 


00 

z(x. ,0)  ■ E [zc(x.)  cosm0  + z8(x.)  sinmQ]  (3b) 

~ i - — m i -m  l 

m-0 

o 

in  which  z(xq,0)  is  the  state  vector  z at  x-Xq  and  z^Xq)  is  the  vector 

of  coefficients  of  the  Fourier  cosine  series  at  x»x  , etc.  Numerical 

o 

integration  produces  a transfer  matrix  so  that 


4<*i> 

4lxi)  ’um4(x0) 


(4) 


c s 

in  which  U and  U are  transfer  matrices  for  coefficients  of  the  Fourier 
m m 

cosine  and  sine  series,  respectively.  Then,  _z(x^,0)  can  be  expressed  as 


l(xr 0)-  I ^z^(x0)  cos  me  + U*  z£(x  ) sin  mG]  (5) 

UFO 

To  derive  the  stiffness  matrix,  a direct  relationship  between  £(Xq,9)  and 
z(x^,0)  is  needed.  If  m*0,  1 and  0"O,  0^,  02  are  Introduced  into 
Eqs.  (3a)  and  (5). 


' i<yon 
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w X0'®2^ 
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Q 

r4(v ; 

r4{  v ■ 

a 

Uf  cos0.j 

U*  sfn01 

► S [T2]  < 

4(V  > 

re 

Uf  cos02 

U*  s1n02 

U(vJ 

wJ 

in  which  I • [I]  is  an  Identity  matrix,  0 - [0]  is  a null  matrix, 
and  I cos0  • [I]  cos0^,  etc.  From  Eqs.  6a  and  6b, 


(6b) 


r *(v°n 

L z(x0,62jJ 


IzOc-j.O) 
Z^(x-|  ,0-| ) > 

z/x^  »02)J 


(7) 


Rearrangement  of  the  state  variables  to  express  the  generalized 

forces  in  terms  of  the  generalized  displacements  gives  the  desired 

global  stiffness  matrix.  If  m • 0,1,2,  , n and  0 ■ 0,  0^,  0^»  , 

0.  , [T. ] and  [T-]  are  square  matrices.  However,  if  other  combinations 
zn  l z 

of  m and  0 are  used,  [T^]  and  [T2 ] will  be  rectangular  matrices.  In 
such  cases,  the  pseudoinverse  can  be  used  to  obtain  the  final  equation 
such  as  Eq.  7. 

The  stiffness  matrix  for  region  III  in  Fig.  1 can  be  obtained  in 
a fashion  similar  to  that  employed  above  for  region  I.  Then,  the 
stiffness  matrices  for  region  I and  III  are  combined  with  the  stiffness 
matrix  for  region  II  which  is  obtained  from  the  finite  element  formula- 
tion. Finally  the  resulting  equation  is  solved  for  the  unknown  displace- 
ment variables. 

Derivation  of  the  transfer  matrix  from  the  stiffness  matrix.  - It  was 
shown  in  Ref.  7 that  a transfer  matrix  in  terms  of  boundary  state  vectors 
at  boundaries  x ■ and  x - x^  in  Fig.  1 could  be  derived  from  given  a 

stiffness  matrix  by  condensing  out  the  intermediate  state  vectors 
between  the  boundaries.  That  la,  the  following  relation  can  be  formed. 

1 -■ »"  1 r.  * • 'v  ' * — ~ ■'  — —■■■ " 

^ -a.  1 _ ■ - ■_  *.  . _ T 9 i -f.  - 1 *-  IT.  ' 1 * 1 1 ■■  ■ “ --1  11  11  ' IT” 
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V 


(8) 


The  next  task  Is  to  express  the  state  vectors  In  terms  of  the  Fourier 
coefficients.  For  this  purpose,  use  the  following  relationships. 


I 

I cos01 
I cos02 
0 


I 

I cos9-| 
I cos02 
0 


Introducing  Eq.  9 into  Eq.  8, 


>wVV' 


(10) 


Then,  the  desired  transfer  matrix  for  the  Fourier  coefficients  is 


-o 

.c 


L](x2) 

— 1 ( x2 ) 
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car1  [a]  [a]  |4<x,) 
*4<xi> 


; cu„]  4<y 


4<V 

s, 


MVA1 
1 
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In  which  [Ujj]  Is  the  transfer  matrix  for  region  II. 


Combine  [Ujj]  with  [Uj]  and  [Ujjj]  which  are  the  transfer  matrices 
for  regions  I and  III  respectively. 


- [Um]  [Ujj]  [Uj] 

! 

» M | 

4<v) 

\ 

Um 

l 

1 4<xo7 

I 

4<V 

1 4(xo>l 

) 

1 J 

. 1 ) 

Upon  Introduction  of  the  boundary  conditions,  Eq.  (12)  can  be  solved 
for  the  remaining  unknown  variables. 
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